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First-principles calculations of the conventional and acoustic surface plasmons (CSPs and ASPs) 
on the (111) surfaces of Cu, Ag, and Au are presented. The effect of s — d interband transitions on 
both types of plasmons is investigated by comparing results from the local density approximation 
and an orbital dependent exchange-correlation (xc) potential that improves the position and width 
of the d bands. The plasmon dispersions calculated with the latter xc-potential agree well with 
electron energy loss spectroscopy (EELS) experiments. For both the CSP and ASP, the same trend 
of Cu<Au<Ag is found for the plasmon energies and is attributed to the reduced screening by 
interband transitions from Cu, to Au and Ag. This trend for the ASP, however, contradicts a 
previous model prediction. While the ASP is seen as a weak feature in the EELS, it can be clearly 
identified in the static and dynamic dielectric band structure. 



PACS numbers: 73.20.Mf, 73.21.-b, 71.45. Gm 

The collective electronic excitations at surfaces, known 
as surface plasmons, have generated extraordinary in- 
terest in the past decade due to the diverse potential 
applications in sensing, imaging, surface enhanced spec- 
troscopy, catalysis and solar energy harvesting^. Noble 
metals such as gold and silver are widely used in exper- 
imental plasmonics due to their stability and well con- 
trolled plasmonic properties. Besides the widely inves- 
tigated 'conventional' surface plasmons (CSPs), which 
have finite plasmon energies for all wave vectors^, the 
noble metal surfaces support another type of surface exci- 
tations, namely the acoustic surface plasmons (ASPs) 3 -^. 
The ASP exhibits a linear dispersion of plasmon energy 
(w aS p) as a function of wave vector with cj asp — > as 
q — > 0. It was predicted to exist on surfaces supporting 
Shockley surface states within the bulk electronic band 
gap and in a two dimensional electron gas on top of a 
substrate^'. Experimentally, the ASP has been identi- 
fied on Be(0001) 5 , Cu(llip and Au( 11 1^ surfaces, as 
well as in graphene adsorbed on substrates^!. 

Theoretically, the CSP and ASP of noble metal sur- 
faces have been treated using two different models. The 
CSP was modeled using a jellium surface, taking the 
screening by d electrons into account via an effective di- 
electric constantPS. For the ASP, a one-dimensional (ID) 
model potential was constructed to reproduce the main 
features of the Shockley surface states (SS) and their un- 
derlying bulk stated. However, these model calculations 
do not account for the presence of d states and are not 
fully ab-initio. 

While experiments are so advanced that not only 
the plasmo ns o n the clean noble metal surfaces can 
be measured^, but their behavior in the presence of 
adsorbates^ and disorder^ has also be studied, a uni- 



fied ab-initio treatment of both kinds of plasmons is yet 
missing. 

In this paper, we use time-dependent density func- 
tional theory with the adiabatic local density approxi- 
mation (ALDA) to calculate the CSP and ASP of the 
(111) surfaces of the noble metals Cu, Au and Ag. 
Single-particle energies and orbitals are obtained using 
the GLLBSC xc-potential to get quantitatively correct d 
band positions and plasmon energies - 3 4 The calculated 
dispersions for both the CSP and ASP are in good agree- 
ment with available experimental EELS measurements 
on Cu and Au surface, despite the fact that the calcu- 
lated ASP energies for Au(lll) are underestimated by 
more than 0.5 eV for q > 0.2 A -1 . While the CSPs are 
strongly influenced by the s-d interband transitions, the 
ASPs are only slightly redshifted. For the ASP energies, 
we found a trend of Cu<Au<Ag, contradicting the trend 
derived from previous model calculations Ag<Cu<Au 4 . 
Finally, we demonstrate that the ASP can be clearly iden- 
tified from the static and dynamic dielectric band struc- 
ture of the noble metal surfaces. 

All calculations were perfor med using the projector- 
augmented wave method GPAV^EEScode. The (111) sur- 
faces were modeled by slabs of either 10 or 24 atoms 
thickness (ca. 2 or 5 nm) , for the study of CSP and ASP, 
respectively. The thicker slabs are essential to converge 
the ASP in the long wave length limit where interaction 
between the plasmons on the two sides of the slab be- 
comes more significant. For comparison, the ground state 
calculati ons w ere performed using both the LDA and 
GLLBSCPESI functionals. The plasmon energies were 
obtained as the peaks in the loss function 

S(q, w) = -Imegio G , =0 (q, w), (1) 
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FIG. 1: (Color online) Band structure (left panel) and den- 
sity of states (DOS, right panel) for the (111) surfaces, repre- 
sented by 24 layer slabs (~ 5nm thickness), of (a) Cu, (b) Au 
and (c) Ag, calculated with LDA (grey lines) and GLLBSC 
(black lines). The bands marked in red represent Shockley 
surface states (SS). Zero energy indicates the Fermi level. 



where G and G' are reciprocal lattice vectors and 
£ gg' (l' w ) * s ^ ne inverse dielectric matrix. The latter 
is obtained using linear response time-dependent den- 
sity functional theory^ in the adiabatic local density ap- 
proximation and with single-particle energies and orbitals 
from either the LDA or GLLBSC ground state calcula- 
tions. More details on the calculations can be found in 
RefP 

Figure [l] shows the band structure (left panel) of the 
24 layer slabs used to represent the (111) surfaces of Cu, 
Ag and Au. The corresponding density of states are pre- 
sented in the right panel. For each surface the GLLBSC 
xc-potential lowers the upper edge of the d bands and 
leads to a narrowing of the d band width, compared to 
LDA. The lowering and narrowing of the d bands with 
respect to the LDA calculations were also observed in 
previous GW calculations on bulk CiP and A^3. The 
calculated GLLBSC upper/lower edge of the d bands are 
at around -2/-6 eV (Cu), -2/-8 eV (Au) and -4/-7 eV 
(Ag). These values are in good agreement with angle 
resolved photoemission (ARPES) measurement on poly- 
crystalline Cu and Au surfaces^ and Ag(lll) films^Sl. 
The broadening of the d bands from Cu(lll) to Au(lll) 
(with an increase of the d band width from 4 to 6 eV), 
and the downshifting of the d bands from Au(lll) to 
Ag(lll) (with the upper edge lowers from -2 to -4 eV) 
indicate reduced d electrons screening from Cu, to Au 



and Ag surface. 

Figure [2] (a)-(c) show the CSP and ASP dispersions 
for the (111) surfaces of Cu, Au and Ag, respectively. 
The results of the present work are indicated by squares 
(LDA) and circles (GLLBSC). The CSP is not visible on 
the Cu(lll) surface [panel (a)] due to strong screening 
(Landau damping) of the plasmon by the s — d inter- 
band transitions. On Au(lll) [panel (b)] the CSP ap- 
pears just above 2.5 eV in GLLBSC calculations, but is 
absent in the LDA calculated loss spectrum. The fact 
that no CSP is predicted by LDA is due to the over- 
estimation of the d band position which leads to over- 
screening of the plasmonaEI. On Ag(lll) [panel (c)] the 
CSP appears for q — at around 3 eV and 4 eV for 
LDA and GLLBSC, respectively, and disperses towards 
higher energy for larger q. Overall, this shows a trend 
of Cu(nonexistent)<Au<Ag for the CSP energies. Com- 
paring to the experimental results (crosses), the CSP en- 
ergies obtained with GLLBSC presents a significant im- 
provement over the LDA. This is clearly a consequence of 
the improved description of the d band positions with the 
GLLBSC. In contrast, the ASP dispersions shown in Fig. 
[2] are much less influenced by the d bands, showing less 
than 0.1 eV difference between the LDA and GLLBSC 
results for all surfaces. The ASP energies exhibit the 
same trend of Cu<Au<Ag [see panel (d)], as found for 
the CSP. 

In general, the s — d interband transitions can have two 
effects on the surface plasmons depending on the relative 
energies of the two types of excitations. If the energies 
of the s — d transitions are comparable to that of the 
plasmon, the latter is strongly redshifted and damped. 
In contrast, if the plasmon energy is significantly lower 
than the s — d transitions, only a slight redshift of the 
plasmon energy is observed^. The first scenario clearly 
applies to the CSP while the second describes the case for 
the ASP. While the trend of Cu(nonexistent)<Au<Ag 
observed for the CSP can be explained by the same trend 
in the d band positions (broadening and downshifting), 
it is not obvious why the ASP energies exhibit the same 
trend. 

Before analyzing further the observed trend in the ASP 
energies, we compare our ab-initio ASP dispersions with 
the model predictions of Ref. HI and EELS experiments. 
For Cu(lll) in Fig. | (a), the GLLBSC presents a 
slight improvement over the model by predicting some- 
what lower ASP energies, although in this case both the 
model and LDA are already in quite good agreement 
with experiment^. These ASP energies lie just above 
the upper boundary of the electron-hole (e-h) pair con- 
tinuum, the latter agrees with previous calculations^. 
For Au(lll), the GLLBSC result for the ASP agrees well 
with experiments at small momentum transfer, but de- 
viates at larger q exhibiting a too small slope for the 
plasmon dispersion. Similar to that for Cu(lll), the 
calculated ASP dispersion for Au(lll) coincides with 
the upper boundary of the Shockley e-h pair continuum. 
Since Au(lll) undergoes a surface reconstruction with a 
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FIG. 2: Conventional and acoustic surface plasmon (CSP and ASP) dispersions for the (111) surfaces of Cu, Au, and Ag 
calculated with LDA (hollow squares) and GLLBSC (solid dots). The plasmon energies are obtained as peaks in the loss 
function 5*(q, oj), the examples of which are shown in Fig. [4] as black lines. Previous model results (dashed line) and EELS 
experimental data (crosses)^ are shown for comparison. Shaded areas mark the Shockley electron-hole pair continuum calculated 
with GLLBSC. (d) Summary of the ASP dispersion obtained with the GLLBSC potential for each of the three surfaces. 



22 x \/3 unit celP, we investigated the effect of strain by 
compressing the Au lattice constant by 4.4%. Insignif- 
icant changes (around 0.1 eV blueshift) were found for 
the ASP energies. A possible reason for the discrepancy 
could be the spin-orbit coupling, which splits the sur- 
face state bands into two subbands in k-space^, but is 
out of scope of this work. For Ag(lll), the GLLBSC 
ASP dispersion has a significantly larger slope than the 
model prediction and is around 0.2-0.4 eV (corresponds 
to q = 0.1-0.2 A" 1 ) above the e-h pair continuum. A 
previous ab-initio calculations on the ASP of a Be(0001) 
surface also shows that the calculated ASP lies slightly 
above the e-h pair continnunP. To the best of our knowl- 
edge, no experimental results are available for the ASP 
on Ag(lll). Compared to Au(lll), the Ag(lll) surface 
does not undergo a surface reconstruction (although eas- 
ily oxidized) and spin-orbit coupling effects are negligible. 
Consequently, experimental results for Ag(lll) would be 
highly important to advance our understanding of the 
surface plasmons on noble metal surfaces. 

The trend of Cu<Au<Ag predicted for the ASP en- 
ergies contradicts the model prediction of Ag<Cu<Au. 
In the model, the ASP energies are directly correlated 
with the Fermi velocity of the Shockley surface band. 
By rigidly shifting the surface bands in energy relative to 
the Fermi level (by up to 0.4 eV) in our ab-initio calcula- 
tions we are able to keep the effective mass of the surface 
bands unchanged while altering the binding energy and 
thus the Fermi velocity of the surface states. However, 
we found that the ASP energies have very weak depen- 
dence on the binding energy and the Fermi velocity of 
the surface states. On the other hand, according to Fig. 
[2j the ASPs seem to be highly correlated with the upper 
edge of the e-h pair continnum. However, the upper edge 



of the e-h pair continuum differ by at most 0.1 eV (at 
q = 0.2 A" 1 ) among the three metal surfaces and thus is 
not sufficient to account for the calculated difference in 
the ASP energies at the same q. Furthermore, the dis- 
crepancy, shown in Fig. [2](c), between the ab-initio and 
the model calculations for the e-h pair continuum of the 
Ag surface, suggests that a ID harmonic approximation 
for the surface bands might not be sufficient in describ- 
ing the e-h pair continuum in the case of Ag. In order to 
understand the trend of the ASP energies, we introduce 
in the following the static and dynamic dielectric band 
structure of the surfaces. 

Fig. [3] (a)-(c) show the static dielectric band struc- 
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FIG. 3: (Color online) Static dielectric band structure of 
the (111) surfaces of Cu, Au, and Ag represented by 24 layer 
slabs (~ 5nm thickness). e~ are the real eigenvalues of the 
static inverse microscopic dielectric matrix e^Q,. The bands 
marked in red represent the surface modes (SMs). 



4 




w (eV) 



FIG. 4: (Color online) Imaginary part of the dynamic dielec- 
tric eigenvalues — Ime^ 1 (q, lj) for the surface states (i = SM, 
red lines in Fig. [3]l plotted as a function of cu at |q| = 
0.125/A for (a) Cu(lll), (b) Au(lll) and (c) Ag(lll), respec- 
tively. They are compared to the macroscopic loss function 
5(q, lj) (black lines) at the same q. 

ture e~ 1 (q) for the 24 layer slabs representing the (111) 
surfaces of Cu, Au and Ag, respectively. Following Ref. 
[26l e~ 1 (q) are obtained as eigenvalues by diagonalizing 
its inverse dielectric matrix £qq, (q, lj — 0) at each q 
in the irreducible Brillouin zone (BZ) . The dielectric 
band structures of all three metal slabs consist of two al- 
most continuous groups of bands separated by a gap from 
around 0.2-0.5 (Cu), 0.3-0.7 (Au) and 0.4-0.6 (Ag). The 
density of the individual bands lying within these groups 
increases for thicker slabs and are consequently related 
to the bulk. In contrast, the single band marked by red 
lying inside the gap of each surface does not change as 
the slab is made thicker. It implies that this band relates 
to the surface state of the slab. To verify, we checked 
the corresponding eigenvector of the red band and found 
it localized at the surface of the slab decaying exponen- 
tially into the slab, which is the characteristic of a sur- 
face state. In fact, the red band is rather two degenerate 
bands in agreement with the fact that the slab contains 
two identical surfaces. Furthermore, we have calculated 
the dielectric band structure of an Al(lll) 20 layer slab 
which docs not support electronic surface states. The 
dielectric band structure (not shown) presents only one 
group of dense bands, i.e. there is no gap, and no bands 
that can be related to the surface. Therefore, the exis- 
tence of the red band (from hereon denoted the surface 
mode (SM)) is a consequence of the SS in the electronic 
band structure. While the electronic SS (red lines in Fig. 
[I]) exist only at small q and merge into the bulk bands 
at larger q, the SMs (red lines in Fig. [3]) are much more 
prominent and are visible inside the bulk dielectric gap 
up to the zone boundary. 



We have argued above that the SM in the dielectric 
band structure is a consequence of the SS in the electronic 
band structure. In the following we present further proof 
of this statement. Figure [4] shows the imaginary part of 
the dynamic dielectric eigenvalues — Ime^~ (q, lj) for the 
surface mode (i — SM) at a particular |q| = 0.125/A. 
The dynamic dielectric eigenvalues are calculated for 
each lj following the same diagonalizing procedure as for 
the static case. The macroscopic loss functions S(q, lj), 
which can be directly compared to the EELS spectra, are 
plotted as well. For each surface, — Imeg M (q, lj) resem- 
bles the loss function and exhibits a peak at the ASP 
energy. Such a peak in — Imer (q, lj) is absent for all 
the other bands (i ^ SM) in the dynamic dielectric band 
structure. This suggests that the ASP is indeed originat- 
ing from the SM. 

Finally, we propose an explanation for the trend of 
Cu<Au<Ag observed for the ASP energies. Shown in 
Fig. [3J for each q, the value of e^ 1 corresponding to 
the SM exhibits a trend of Cu<Au<Ag. Recalling that 
e^ 1 gives the total potential due to an externally ap- 
plied potential, this shows that the electronic screening 
for the SM (thus ASP) is strongest on Cu(lll), weaker 
on Au(lll) and weakest on Ag(lll). Again, this is a 
consequence of the position and width of the d bands 
as discussed in connection with the CSP. Furthermore, 
the existence of the SM within the bulk dielectric gap up 
to the BZ boundary agrees with ASP being observed at 
relatively large q. 

In conclusion, the conventional and acoustic surface 
plasmons on the (111) surfaces of Cu, Au and Ag were 
investigated using time-dependent density functional the- 
ory with the adiabatic local density approximation. Our 
ab-initio results agree well with available EELS experi- 
ments. The energies of both conventional and acoustic 
surface plasmons follow the same trend of Cu < Au <Ag. 
The trend is attributed the reduced d bands screening 
from Cu, to Au and Ag. Experiments of the ASP en- 
ergies on Ag(lll) surface will be very important to test 
our findings and further advance our understanding of 
plasmons and electronic screening at metal surfaces. 
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